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A FURTHER COMPARISON OF 50-MILLIBAR GEOPOTENTIALS 
OBTAINED FROM SATELLITE INFRA-RED SPECTROMETER 
SOUNDINGS AND ANALYSED CHARTS 


By V. BAILEY 


SUMMARY 


In March 1975 the National Oceanic and Atmospheric Administration (NOAA) introduced 
a ‘regression’ method of retrieving Satellite Infra-red Spectrometer (SIRS) temperature 
profiles from Vertical Temperature Profile Radiometer (VTPR) data. Comparisons of 50 mb 
SIRS geopotentials with those derived from a series of 50 mb charts, based upon radiosonde 
data, have been carried out for the period 25 August 1975 to 12 April 1976. Although the 
quality of the SIRS data is still rather variable, there has been a marked reduction in bias 
errors and some decrease in the standard deviation of the SIRS minus analysed differences 
since the new retrieval method became operational in August 1975, especially over the 
Atlantic and southern North Pacific Oceans. However, north of 50°N in the Pacific, a bias 
error persists. 


In a previous comparison of 50 mb geopotentials from Satellite Infra-red 
Spectrometer (SIRS) soundings and analysed charts, which covered the period 
June 1974 to May 1975 (Watson and Bailey, 1976), the quality of the SIRS data 
was found to be very variable. The mean SIRS minus analysed differences 
displayed a marked latitude dependence, the SIRS geopotentials being com- 
parable to the analysed geopotentials at 60°N and about 10 geopotential 
decametres too high between 20° and 30°N. The mean difference varied with 
time, with discontinuities whenever the instrument or spacecraft, or both, 
changed. It was noted in that paper that a new technique to produce the SIRS 
thicknesses from observed Vertical Temperature Profile Radiometer (VTPR) 
radiances was introduced by the National Oceanic and Atmospheric Adminis- 
tration (NOAA) on 13 March 1975. In this technique (Hayden, 1976) each 
thickness is derived from a regression equation involving the deduced clear- 
column radiances for several VTPR channels. 

The regression coefficients were derived from an analysis of co-located 
radiosonde and SIRS observations. These coefficients were initially calculated: 
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from comparisons made during January and February 1975 and subsequently 
updated on 21 August 1975 with collocations made in June and August. It is 
understood that, since this date, there has been no change of coefficients or 
methods. Instrument No. 2 on NOAA 4 was in use throughout the period of this 
study. Indeed, it has been learnt that, contrary to Figure 1 of Watson and 
Bailey (1976), instrument No. 2 on NOAA 4 was brought into use from Decem- 
ber 1974. Thus the abrupt decrease in mean differences in March 1975 was 
solely the consequence of introducing the regression method for deducing SIRS 
thicknesses. 

For this present investigation, comparisons have been made on 34 occasions, 
seven days apart, during the period 25 August 1975 to 12 April 1976. The 
method of comparison was similar to that used in the previous investigation; 
values were taken from a daily (00 GMT) series of northern hemisphere 50 mb 
charts, drawn up using radiosonde data only. Comparisons were made for SIRS 
observations made within -+-3 hours of 00 GMT and therefore were confined to 
the Atlantic and western Pacific Oceans. The average number of SIRS minus 
analysed differences on each occasion was 33, with a maximum of 47 and a 
minimum of 15. It was found that the two populations of differences for the 
Atlantic and Pacific respectively were significantly different at the 1 per cent 
level. Figure 1 shows the mean differences, together with the standard errors of 
the means, at intervals of seven days for both areas. The average of the mean 
differences for the Atlantic was + 1-0 geopotential decametres and for the Pacific 


ATLANTIC 
Average mean difference 1-0, 
SD of mean differences 35 
geopotential decametres 


PACIFIC 
Average mean difference 6:5, 
SD of mean differences 43 
geopotential decametres 


Differences (geopotential decametres) 











FIGURE 1—MEAN (SIRS MINUS ANALYSED) GEOPOTENTIAL DIFFERENCES AT 
50 MILLIBARS FOR 25 AUGUST 1975 To 12 APRIL 1976 AT INTERVALS OF 
SEVEN DAYS FOR ALL AVAILABLE COMPARISONS BETWEEN 20° AND 70°N 


Vertical bars indicate standard errors of means. 
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-+-6°5 geopotential decametres. The increased variability of the mean differences 
for the winter months when compared with autumn is particularly noticeable. 
For the Atlantic area, an examination of the data for the period August 1975 
to April 1976 in 10° latitude bands showed (bottom plot of Figure 2) the varia- 
tion of the mean difference with latitude to be small. It varied from 1-8 geo- 
potential decametres between 20° and 30°N to —0-6 between 40° and 50°N and 
2:8 between 60° and 70°N. In order to assess the impact of the change of 
retrieval method directly, data for part of the previous study period (15 December 
1974 to 13 March 1975) have been reanalysed. During the latter period, the same 
instrument (No. 2) on NOAA 4 was being used but SIRS values were being 
derived by the ‘minimum information method’. These results are included in 
Figure 2. The current regression method of retrieval has considerably suppressed 
the variation with latitude and eliminated the systematic error at the 50 mb level. 


e— ——e Atlantic, old (minimum information) method 
° e Atlantic, new (regression) method 
x — =x Pacific, old method 
x—— Pacific, new method 
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FIGURE 2—VARIATION WITH LATITUDE OF THE MEAN DIFFERENCES OF 
GEOPOTENTIAL BETWEEN SIRS AND ANALYSED DATA 


Vertical bars indicate standard errors of means. 
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The data for these two periods are summarized in Table I. N is the number of 
observations. 


TABLE I—50-MILLIBAR SIRS MINUS ANALYSED GEOPOTENTIAL DIFFERENCES 
OVER THE ATLANTIC (IN GEOPOTENTIAL DECAMETRES) 


Minimum information method Regression method 
15/12/74-13/3/75 25/8/75-12/4/76 

Mean SD N M 

60-70°N 9-4 
50-60°N 11-4 1 
40-50°N 13-3 i 
30-40°N 17-3 1 
20-30°N 20-6 1 


4 104 
0 203 
9 239 
9 339 
3 260 


PAWMHOYV 
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9. 
0: 
0- 
3 
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Using the standard deviation (SD) as a measure of variability of the differences, 
the Table suggests that there has been a considerable improvement in the quality 
of the SIRS data at 50 mb over the Atlantic at lower latitudes. 

For the Pacific (Table II) there is still a mean systematic error. The variation 
with latitude is also more marked, the SIRS minus analysed difference increasing 
from 4-1 geopotential decametres between 20° and 30°N to 11-3 geopotential 
decametres between 50° and 60°N. (Only 14 observations were available north 
of 60°N.) However, the latitude dependence (Figure 2) shows a considerable 
improvement south of 50°N. The standard deviations are also improved in this 
region and are only slightly larger than for the Atlantic. 


TABLE 1I—50-MILLIBAR SIRS MINUS ANALYSED GEOPOTENTIAL DIFFERENCES 
OVER THE PACIFIC (IN GEOPOTENTIAL DECAMETRES) 
Minimum information method Regression method 


15/12/74-13/3/75 25/8/75-12/4/76 
Mean D N Mean 


60-70°N — 14 16°6 
50-60°N ; : 166 11:3 
40-S0°N : : 193 : 

30-40°N 15- F 227 
20-30°N 23° : 267 


The bias in the ~ -th of the Pacific area still persists, but it may be asked 
whether this reflects an error in the SIRS values or in the analysis. The method of 
analysis requires that corrections be applied to the various types of sonde 
observations to make them compatible with each other, using the Kew Mk 2b 
radiosonde as a standard. These corrections, obtained by the comparisons of 
sondes with the objective analysis, are small (up to —4 geopotential decametres) 
over the Pacific in low latitudes but fairly large (up to —11) in high latitudes for 
00 GMT (when this area is in sunlight). Campbell and May (1976) have carried 
out an investigation using co-located SIRS/sondes, analysing their results by 
types of sonde. In relation to the Kew sonde, they found differences which, 
while having the same sign as the corrections that we apply in our chart analysis, 
were smaller. For instance, they found that the mean difference (in geopotential 
decametres) between the Kew sonde and the Russian sonde was approximately 3 
for the 1000-100 mb thickness, which would correspond to about 4 for the 
1000-50 mb thickness, compared with the correction of 11 applied during chart 
analysis. Use of the Campbell and May figures would greatly reduce the latitude 
dependence in the north of the Pacific area. This, of course, assumes that SIRS 
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values are in some way more consistent than radiosonde values in judging the 
difference between sondes at night over the Atlantic and in daylight over the 
Pacific. The only sonde for which Campbell and May were able to investigate 
differences between day and night was the Japanese sonde. Their findings were 
that the SIRS minus analysed values were 3-9 geopotential decametres greater 
by night than by day, in close agreement with 4-7 obtained from the objective 
analysis scheme. This result suggests that SIRS soundings for a given location do 
not show any significant day-to-night difference. However, other possibilities 
exist in the SIRS analysis for the introduction of bias between the Atlantic and 
the Pacific. NOAA uses several sets of regression coefficients, selected by values 
of observed clear-column radiances in the window channel and a channel (No. 2) 
with a weighting function in the stratosphere. The climatic difference between 
the North Pacific and North Atlantic temperatures at low stratospheric levels 
(Labitzke et alii, 1972) presumably causes a bias towards different sets of 
coefficients for the production of SIRS values for the North Pacific and North 
Atlantic. This is something that the ordinary user of SIRS soundings cannot 
check without also accessing the clear-column radiance data. 

In conclusion, it is clear that there has been a marked improvement in quality 
of SIRS geopotential estimates at 50 mb since the regression method of retrieval 
was introduced, mainly because of the reduction in bias error for the Atlantic 
and southern North Pacific. For the northern part of the Pacific area our analysis 
suggests either that there is a bias error in some of NOAA’s sets of regression 
coefficients or that our present estimates of sonde corrections over the north of 
the Pacific in daylight are misleading. Continuation of the present form of 
comparison of SIRS values with the analysed charts is now no longer valid. 
This is because the 50 mb charts are constructed from a 100 mb chart, produced 
in the Central Forecasting Office, which now contains an element of SIRS 
observations. Hooper (1975) has suggested that satellite radiances calculated 
from upper-air soundings may provide a useful basis for comparison of sondes. 
Studies of SIRS minus analysed differences will be continued in order to improve 
the understanding of both SIRS and sonde behaviour. 
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A NUMERICAL INDEX TO MONITOR THE AFRO-ASIAN 
MONSOON DURING THE NORTHERN SUMMERS 
By J. FINDLATER* 


SUMMARY 


The identification of the core of the major low-level air current of the Afro-Asian summer 
monsoon at a topographically fixed position over eastern Africa has allowed an attempt to be 
made to monitor the low-level flow and relate it to rainfall downstream over western India. 
An index of the southerly flow over eastern Africa at a station a few degrees south of the 
Equator for the month of July, for a period of 24 years, is compared with the July rainfall of 
ten stations in the western part of the State of Maharashtra, India. 

It is found that months with an index of high or low wind correspond well with months of 
high or low rainfall respectively, especially when two-year overlapping averages are used. 
An interesting feature of the analysis is that there is a one-year lag between maxima and 
minima of the wind index over eastern Africa and the corresponding features of the rainfall of 
western India. A tentative calculation is made to illustrate how the lag might be used for 
long-range rainfall prediction. 


1. INTRODUCTION 


It has been reported earlier (Findlater, 1969a, 1969b and 1971) that the low-level 
airflow of the Afro-Asian summer monsoon is organized into a narrow high- 
speed current circulating at about 1-5 km above mean sea level in the western 
periphery of the monsoon system. A simplified map of the mean flow at 1 km 
for the month of July is shown in Figure 1, and it is noticeable that the core of the 
current, though primarily an oceanic phenomenon, passes over the flat arid lands 
of eastern Kenya, eastern Ethiopia and Somalia. 

It was over these land areas that low-level jet streams, with speeds reaching 
25-50 m/s, were first noticed some years ago (Findlater, 1966 and 1967). These 
low-level jet streams, or surges in the major current, have been located over the 
ocean as well as over land—usually in the vicinity of the axis of maximum flow 
shown in Figure 1 and orientated along it. 

The fact that the axis of the major current passes over land was used in an 
earlier analysis to monitor the flow at an equatorial station in eastern Africa and, 
when five-day overlapping means were used, a close relationship between 
pulsations of the flow over eastern Africa and fluctuations in the rainfall of 
western India was noticeable (Findlater, 1969a). Similar but unpublished 
analyses for some other years have confirmed this relationship. It should be 
noted, however, that the correspondence between cross-equatorial flow over 
eastern Africa and the rainfall of western India can only be located if the upper- 
wind stations which are used lie under or very close to the core of the current 
shown in Figure 1. Attempts to use radar-wind data from Nairobi (01°18’S, 
36°45’E, 1798 m above mean sea level) or Dar es Salaam (06°53’S, 39°12’E, 
55 m) have proved fruitless because Nairobi is not affected by the current and 
Dar es Salaam lies on the fringe of it and is distant from the core. Recourse 
must be had to pilot-balloon data, and the earlier analysis used data from 
Garissa (00°29’S, 39°38’E, 128 m), a station almost on the Equator and close to 
the core of the current. 





* At present on secondment from the Meteorological Office to the East African Meteoro- 
logical Department, Nairobi, Kenya. 
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FIGURE 1—MEAN MONTHLY AIRFLOW AT 3000 ft (1 km) IN JULY 
M—Mombasa G—Garissa 


(The dots over India represent the ten selected rainfall stations.) 


The aim of this present study is to monitor the flow over eastern Africa, using 
mean values for the month of July from as many years as possible, to compare 
the vigour of the monsoon from year to year and to relate it to the mean July 
rainfall downstream over western India. 


2. SELECTION OF A MONITORING PARAMETER 


Although the pilot-balloon station at Garissa is ideally situated for monitoring 
the flow at the core of the major current where it crosses the Equator, upper- 
wind soundings commenced there only in 1962. However, another station— 
Mombasa (04°02’S, 39°37’E, 57 m)—lies close to the position where the core of 





172 Meteorological Magazine, 106, 1977 


the current first comes inland from the Indian Ocean, and the length of record of 
upper winds is much greater. Also, the sounding program at Mombasa has been 
much more regular than at Garissa. For these reasons Mombasa was selected 
as the monitor station for this study. The mean characteristics of the low-level 
flow at these two stations is shown in Figure 2 where the mean southerly and 
easterly components of the July winds, up to the 3 km level, illustrate the jet-like 
profile of the southerly wind. The peak value at Garissa is higher because that 
station lies a little nearer to the core than Mombasa. 
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E’ly components S'ly components "~~ ~* Garissa (1962-70) 
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FIGURE 2—MEAN COMPONENTS OF THE LOW-LEVEL CURRENT IN JULY AT 
MOMBASA AND GARISSA 


The index of cross-equatorial flow is the mean of the monthly values of the southerly 
component at 1000, 2000, 3000, 4000 and 5000 ft above MSL. 
For Mombasa the 22-year mean of the index is 20-0 kn. 


Points A and B represent the components of the mean daily surface wind in July at 
Mombasa, as calculated from data given by Ramsey (1971). 


When monitoring the flow each July only the lower half of the profile shown in 
Figure 2 is used because many more pilot balloons reach the 1-5 km level than 
the 3 km level. The monitor value is calculated for each July as the mean of the 
southerly components of the wind at the levels of 0-3, 0-6, 0-9, 1-2 and 1-5 km 
above MSL (corresponding to winds measured at 1000, 2000, 3000, 4000 and 
5000 ft above MSL). The 22-year mean of the monitor value for Mombasa for 
the years 1953 to 1974 is 20-0 kn (=10-3 m/s) and it is shown graphically in 
Figure 2. The monitor values for July in each of the years 1953-76 are listed in 
Appendix I*. Extreme values of 16-9 and 22-6 kn have been recorded. The 
monitor value for each year, hereinafter referred to as the wind index, is also 
shown in (a) of Figure 3. 





*Pilot-balloon winds have been measured at Mombasa since 1937 but all daily records 
prior to 1953 have been destroyed. Only an incomplete set of monthly frequency tables for a 
ew upper levels remains. 
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FIGURE 3—CROSS-EQUATORIAL FLOW AT MOMBASA IN JULY (a), AND TOTAL 
JULY RAINFALL AT THREE STATIONS IN WESTERN MAHARASHTRA (b), (c) AND (d) 


(b) represents a high-rainfall station, (c) a medium-rainfall station, and (d) a low-rainfall 
station. 


3. COMPARISON OF LOW-LEVEL FLOW OVER EASTERN KENYA AND THE 
RAINFALL OF WESTERN INDIA 


The July rainfall of western India varies considerably from year to year and (b), 
(c) and (d) of Figure 3 show the total July rainfall at three stations in the western 
part of the State of Maharashtra, India for the years 1953-76. The three stations 
represent a high-rainfall station (Bombay), a medium-rainfall station (Kolhapur) 
and a low-rainfall station (Poona). The 1953-76 mean July rainfall values at 
these stations are 873-0, 390-0 and 208-2 mm respectively. Although these 
stations are at different altitudes and have different exposures, it is clear that all 
lie in the same rainfall regime. Wet and dry years are generally reflected at all 
three stations. 

Data from these three stations can be combined for comparison with the 
Mombasa wind index. This comparison is shown in Figure 4. The profiles of 
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FIGURE 4—COMPARISON OF THE INDEX OF CROSS-EQUATORIAL FLOW OVER 

EASTERN AFRICA IN JULY (a), THE MEAN JULY RAINFALL AT THREE STATIONS 

IN WESTERN INDIA (b), AND THE MEAN JULY RAINFALL AT TEN STATIONS IN 
WESTERN INDIA (c) 


The dots represent the mean July value for two consecutive years. The crosses in (a) represent 
the two-year mean July values at Garissa, for comparison with those of Mombasa. 


the wind index, (a) of Figure 4, and of the mean rainfall at the three stations in 
India, (b) of Figure 4, show some similariues and also some differences. To 
assist interpretation, two-year overlapping mean values have been calculated. 
These values, plotted as dots joined by broken lines; illustrate a close association 
between the two derived parameters. Generally, peaks and troughs are nearly 
coincident but the striking feature is the lag between the two curves of overlapped 
mean values. Peaks and troughs in the rainfall lag behind, by one year, the 
corresponding features of the wind-index curve. 

The physical causes of the lag are not yet known but it is hypothesized that the 
strength of the wind index influences the ocean currents and surface tempera- 
tures of the Arabian Sea in such a way that sea-surface temperature anomalies 
persist for one year, thus affecting the stability of the air, evaporation from the 
sea surface, and hence the rain-bearing potential of the Indian summer monsoon 
air. This hypothesis requires verification. 
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In (a) of Figure 4 crosses have been inserted to represent the two-year over- 
lapping means of the wind index at Garissa in order to illustrate the corres- 
pondence between upper-wind stations near the core of the low-level current. 
The averages for Garissa have not been calculated for the years after 1968-69 
because of paucity of data. For example, in 1970 upper-wind measurements at 
Garissa were made on only nine days in July. Nevertheless the Garissa wind 
index for earlier years shows the generally higher speeds recorded at that station, 
and the curve of the two-year averages shows the same trends as at Mombasa. 

The comparison of the Mombasa wind index and the Indian rainfall is 
extended by using a total of ten rainfall stations in western Maharashtra instead 
of only three; these data are shown in (c) of Figure 4. Again the same general 
correspondence of curves is evident, though less pronounced, and the one-year 
lag is seen in both smoothed and unsmoothed values. 

All wind and rainfall data used in these analyses are tabulated in Appendix I. 
Station details are listed in Appendix IV. 

The correspondence between the rainfall curves of (b) and (c) of Figure 4 
is not solely due to the three rainfall stations of (b) being included in the ten 
stations used to construct (c), since all stations used in this study show similar 
features in their rainfall regimes. For example, if another three stations, Dahanu, 
Harnai and Ratnagiri, are analysed in similar fashion, the pattern is the same and 
the lag of one year between troughs and ridges in the wind index and rainfall 
curves is repeated. It should be noted, however, that the relationships reported 
in this paper do not necessarily apply to rainfall stations which lie outside the 
area selected for study, or to other months. 


4. RAINFALL PREDICTION 


The lag of one year between significant features of the wind and rainfall, shown 
in Figure 4, suggests that the curve of Indian rainfall may be extrapolated one 
year ahead for the month of July only. The trend or change of the two-year 
mean values of the wind index may be determined and related to the trend of the 
two-year mean rainfall of western India one year later; for example, the change 
in the two-year mean wind index between 1953-54 and 1954-55 is —1-6 kn, and 
is related to the change in the two-year mean rainfall for the ten stations in 
western India for 1954-55 and 1955-56, —79-2 mm. These comparative values, 
listed in Appendix III and plotted in Figure 5, are calculated from data given in 
Appendix I. The line of best fit in Figure 5 is inserted by the method of least 
squares, whilst parallel and equidistant lines contain the range of observations. 
The correlation coefficient between the two-year wind-index change and the two- 
year mean rainfall change one year later is +0-82. Values for the correlation 
coefficient (r) for a lag varying from —2 to +2 years are: 
Lag (years) —2 —1 0 +1 +2 
r —0-43 —0:54 +0-41 +032 +0-04. 

Provided that the trends of wind and rainfall over the last 24 years are pre- 
served in future years, and provided that the pilot-balloon data at Mombasa 
remain sufficiently regular to furnish an accurate measure of the vigour of the 
cross-equatorial current (and bearing in mind the inherent uncertainty of 
extrapolating time-series into the future), then the type of diagram shown in 
Figure 5 might be used for experimental predictions of the July rainfall of 
western Maharashtra for one year ahead. 
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For example, the 1974-75 and 1975-76 mean values of the Mombasa wind 
index are 18-7 and 18-8 kn respectively, a change of +0-1 kn. From Figure 5a 
change of +-0-1 kn corresponds to a change of +-2-0 mm (+80 mm) between the 
1975-76 and 1976-77 mean July rainfall for the ten selected stations in western 
India. Because the 1975-76 mean rainfall is known to be 535-4 mm, the 1976— 
77 mean value may be predicted to be 535-4 + 2:0 = 537-4 mm (+80 mm). 

Furthermore, the 1976 measured value of the rainfall is 551-1 mm and the 
predicted mean for 1976-77 is 537-4mm, therefore the predicted mean for 
July 1977 is: 

537-4—(551-1—537-4) = 523-7 mm (+160 mm). 


This predicted amount represents 97 per cent (29-6 per cent) of the 24-year 
average July rainfall (540-1 mm) for the ten stations in western Maharashtra. 
The calculation of the forecast value may also be made graphically by using 
Figure 5 and (c) of Figure 4. 


8 


Rainfall excess (mm) 





r=082 
y=4143x—212 


~ 
| 
x 
2 
me] 
£ 
me) 
P 
> 
mo] 
= 
= 
vo 
ae) 
o 
= 
. 
ie] 
o 
~ 
Pe. 
<= 
= 
= 
2 
me] 
= 
c 
- 
2 
“” 
o 
2 


& 
ww 
c 
© 

m4 

= 
o 

2 
“ 

o 

= 

- 
o 

<) 

oe 
= 
re) 
£ 
c 
5 
o 
Ee 
. 
5 
oO 
> 
f 
nN 
= 
og 
D> 
Cc 
re) 
xo 
UO 


Rainfall deficit (mm) 
| 
S 
o 











| | l | | 
—25 —20 -15 -—10 —05 0 0-5 10 15 20 25 
Change in 2-year mean wind index at Mombesa (kn) —x 





FIGURE 5—CHANGE IN TWO-YEAR MEAN CROSS-EQUATORIAL FLOW OVER 
EASTERN AFRICA RELATED TO THE CHANGE IN TWO-YEAR MEAN RAINFALL 
AT TEN STATIONS IN WESTERN INDIA ONE YEAR LATER 


If similar calculations are made for all years for which data are available, as in 
Appendix II, it is evident that the mean July rainfall of the ten selected stations 
might have been forecast one year ahead with an average error of only 11-3 per 
cent, had the relationships between the Kenya wind index and the Indian rain- 
fall then been known. The maximum error in the forecasts would have been 
29-6 per cent. All details are listed in Appendix II. 

Experimental forecasts using this method were made in August of 1972, 1973, 
1974 and 1975 for the July rainfall of 1973, 1974, 1975 and 1976 respectively, 
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but data from these years have now been incorporated into the dependent data 
set and the correlation coefficient and line of best fit have been recalculated to 
the values now shown. However, it will be many years before the forecasting 
potential of this method can be adequately tested against an independent data 
set. 

It is not the intention here to propose an operational forecasting technique as 
such, but rather to demonstrate that careful monitoring of the flow near the core 
of the cross-equatorial current over eastern Africa and the rainfall downstream 
over India may yield, after much further study, some relationships which might 
be useful in forecasting the rainfall of the Indian summer monsoon. This paper 
is only a short step in that direction. 


5. CONCLUSIONS 


Earlier work relating cross-equatorial flow at low levels over Kenya to the rain- 
fall of western India, using five-day overlapping mean values during the northern 
summer, indicated that surges and lulls in the cross-equatorial flow were reflected 
in the rainfall of western India a few days later. This present study attempts to 
ascribe an index to the core of the cross-equatorial current in July for as many 
years as possible, so that the vigour of the monsoon current may be compared 
between years. 

This index has been found to vary as the mean July rainfall of western 
Maharashtra, especially when two-year averages are used. 

A feature of considerable interest is that there is a lag of one year in maxima 
and minima of rainfall behind those of the wind index. This lag may have some 
long-range predictive value. 
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APPENDIX III—CHANGE IN 2-YEAR MEAN WIND INDEX RELATED TO 
CHANGE IN 2-YEAR MEAN RAINFALL ONE YEAR LATER 


Change in 2-year mean wind 
index at Mombasa, Kenya 


1953/54—1954/55 
1954/55—1955/56 
1955/56—1956/57 
1956/57—1957/58 
1957/S8—1958/59 
1958/59—1959/60 
1959/60—1960/61 
1960/61—1961/62 
1961/62—1962/63 
1962/63—1963/64 
1963/64—1964/65 
1964/65—1965/66 
1965/66—1966/67 
1966/67—1967/68 
1967/68—1968/69 
1968/69—1969/70 
1969/70—1970/71 
1970/71—1971/72 
1971/72—1972/73 
1972/73—1973/74 
1973/74—1974/75 
1974/75—1975/76 


knots 
—1-6 
+0-6 
+0-8 
—0°1 


PUNMeNOoCoon es SS SS 
-ONODAUNW 


Change in 2-year mean rainfall 
at ten selected stations in 
western India, one year later 


1954/55—1955/56 
1955/56—1956/57 
1956/57—1957/58 
1957/58—1958/59 
1958/59—1959/60 
1959/60—1960/61 
1960/61—1961/62 
1961/62—1962/63 
1962/63—1963/64 
1963/64—1964/65 
1964/65—1965/66 
1965/66—1966/67 
1966/67—1967/68 
1967/68—1968/69 
1968/69—1969/70 
1969/70—1970/71 
1970/71—1971/72 
1971/72—1972/73 
1972/73—1973/74 
1973/74—1974/75 
1974/75—1975/76 
1975/76—1976/77 


mm 
—79-2 
+53-6 
—29-9 

onoihey 
—47°5 
+38-1 
+21°2 
—40:1 
—30°8 
+36°5 
+70°5 

—9-3 
—51°5 

—63 
—41-5 

—165-0 
—35-4 

+1120 

+132:3 

+16:3 

—169 
? 


Correlation coefficient (r) = +0°82 (significant at the 1 per cent level of probability). 


APPENDIX IV—LIST OF STATIONS USED IN THE ANALYSIS 


For wind analysis 


Station No. 


63723 Garissa 


Station 


63820 


For rainfall analysis 


43009 
43069 
43001 
43109 
43157 
43192 
43110 
43117 
43003 
43063 


Mombasa 


Ahmednagar 
Baramati 

Dahanu 

Harnai 

Kolhapur 

Panjim 

Ratnagiri 

Sholapur 

Bombay (Santa Cruz) 
Poona 


Position 


00°29’S 39°38’E 
04°02’S 39°37’E 


19°05’N 74°48’E 
18°09’N 74°35’E 
19°58’N 72°43’E 
17°49’N 73°06’E 
16°42’N 74°14’E 
15°29’N 73°49’E 
16°59’N 73°20’E 
17°40’N 75°54’E 
19°07’N 72°51’E 
18°32’N 73°51’E 


Height 
above MSL 
metres 


128 
57 
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551.593.653 


NOCTILUCENT CLOUDS OVER WESTERN EUROPE AND 
THE ATLANTIC DURING 1976 


By D. H. McINTOSH and MARY HALLISSEY 
(Department of Meteorology, University of Edinburgh) 


Table I summarizes the observations of noctilucent clouds (NLC) made over 
western Europe and the Atlantic during 1976 and reported to the Department of 
Meteorology, Edinburgh University. 

Observers’ reports, positive or otherwise, were requested for the months May 
to August. As in previous years, sightings were confined to the period from late 
May to early August. The periods of time during which the clouds were observed 
appear in the second column of the Table. These should not be taken as being 
necessarily the total duration of the display; this is stated where possible, but it is 
obviously difficult, particularly for voluntary observers, to record a display to the 
point of disappearance. Brief notes on the displays appear in the third column. 
In the remaining columns, details of the relevant station co-ordinates are listed 
to the nearest half degree, and the maximum elevation and limiting azimuths of 
the observed cloud, where known. 

Positive reports were received from some 30 stations of the Meteorological 
Office station network of Great Britain and of the Meteorological Service of 
Ireland, ranging from Lerwick to Mount Batten, Plymouth, and from Birr to 
Wattisham; also from a Fair Isle lighthouse-keeper and from 10 voluntary 
observers scattered throughout the UK, Denmark and Norway. We are grateful 
to Drs Gadsden and Jenkins of Aberdeen University for observations recorded 
in their log in connection with their work for the International Magnetospheric 
Study on the polarization of NLC. 

Records of tropospheric cloud amounts were also received from Meteoro- 
logical Office stations, including nights on which no positive observations of 
NLC were made. The confirmed absence of NLC, especially during what is 
statistically their peak incidence time (mid-June to mid-July) is now regarded as 
significant ‘negative’ information. Unfortunately, the incidence of tropospheric 
cloud and haze was rather high during this four week period and no negative 
finding can be regarded as definitely confirmed. 

Time-lapse photography of the clouds was carried out on most nights from 
Edinburgh and revealed one occasion, in particular, on which there was a 
dramatic invasion of the sky, by movement of the clouds from the north, late 
in the night. 

A special study of the display of 18-19 June was made by an observer, 
Dr D. A. R. Simmons of Milngavie, near Glasgow, whose paper on the event has 
been accepted for publication in Weather. 


TABLE I—DISPLAYS OF NOCTILUCENT CLOUDS OVER WESTERN EUROPE AND THE 
ATLANTIC DURING 1976 


Date— Station Max. Limiting 
night of Times Notes position Time elev. azimuths 
UT UT degrees 
27/28 May 0130 Patch of NLC seen from train when 54:-5N 01:5°W 0130 5 
observer travelling NE England 
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TABLE I—continued 


Date— Station 5 Max. Limiting 
night of Times Notes position* elev. azimuths 
UT degrees 


2/3 June 2210-0200 Some doubt by observers in Aberdeen 57°N 02°W 

as to authenticity of NLC to high 

elevation before 2300 but NLC seen 

clearly from more southerly stations 55-5°N 01-:5°W 
one of which later reported veil and 

billow formation. From northerly 55°N 04:5°W 
station NLC faint and diffuse without 54°N 04:5°W 
marked regular structure. At 0020 as 

seen from Aberdeen northern edge of 

cloud field clear of N horizon. 


0115-0150 NLC already at max. when first seen 
at 0115: compact herring-bone- 
pattern patch with N edge 11° NNE- 
NE. Intensity, never high, faded pro- 
gressively until dawn, around 02 


0001-0030 Fairly bright but featureless bands 56°N 03°W 
seen from southern position; later 51°N 02°W 
seen in breaks in tropospheric cloud 
from Edinburgh. 


2110-0215 Bright NLC with fine herring-bone 56°N 03°W 
structure seen from Denmark (photo- 55-5°N 12°E 
graphs available): marked longitu- 55°N 14-S°E 
dinal extension to east; from Liverpool 
seen as veil background to brighter 
bands, fading into haze. 53-5°N 03°W 


2200, 2230 NLC visible through breaks in tropo- 56°N 03°W 
spheric clouds in Denmark; herring- 55-5°N 12°E 
bone pattern discernible Edinburgh. 


2400, 0200 Details obscured by cloud from 54°N 04:5°W 
Ronaldsway; banded structure dis- 53°N 01-5°W 
cernible from Watnall. 


2330 Moderately bright NLC above cloud 56°N 03°W 
bank to 5°. Fairly extensive cloud 
field, but E limit hidden by tropos- 
pheric clouds. 


17/18 June 0045-0200+ Very bright bands and billows with 
veil background became visible after 
midnight from Newcastle. Earlier 
in Denmark, NLC recognized in 
tropospheric cloud breaks. 


Earliest sighting of bright bands from 57:5°N 03-5°W 
Ronaldsway, 1.0.M. Max. extension 57°N 02°W 
southwards around 0200 when 
southern edge of cloud field, though 
obscured by tropospheric clouds at 56:5°N 03°W 
many stations, was overhead at 56°N 04:5°W 
Boulmer, Northumbria, and to eleva- 
tion of 160° from Aberdeen. NLC 56°N 03°W 
structure described variously as 55-5°N 01:5°W 
‘tangled’, ‘feathery’. Measurements 
made from series of photographs 54-5°N 06°W 
taken near Glasgow. 

54°N 04-5°W 


2330-2350 Short-lived period of visibility before | 55-5°N 05°W 
NLC obscured by tropospheric cloud. 


2300 In UK brief recognition through tro- 57°N 02°W 
pospheric cloud break—no identifica- 56°N 03°W 
tion of form. Bright display seen 55-5°N 12°E 
a above tropospheric cloud to 53°N 01-5°W 


0030 In Edinburgh NLC glow visible 56°N 03°W 
pps on tropospheric cloud covering 
N sky to zenith. 
0518-0553 NLC visible from aircraft over New- 48°N 55°W 
foundland. 


23/24 2315 NLC visible through breaks in tro- 56°N 03°W 
pospheric cloud. 





* To nearest 0-5 degree. 
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TABLE I—continued 


Date— : Station Max. Limiting 
night of Times Notes position i elev. azimuths 
UT degrees 


30fJune/1 July 2300-0220 


Faint NLC above cloud bank to 10°. 


Veil with band formation visible from 
all reporting stations (including the 
top of Ben Lomond). Extensive sight- 
ing and detailed report from Alder- 
grove; there, as earlier in Denmark, 
bands, billows and whirl formations 
identified, the tenuous structure 
thickening into broad bands at 
southern edge before fading in bright- 
ening sky. 4° elevation of northern 
edge as seen from Aldergrove. 


Steady, slow-moving display seen 
Ireland and 1.0.M., mainly in NW 
sky, and in tropospheric cloud breaks 
Edinburgh. Weak veil and brighter 
bands, with whirl formation suspected 
around 2300. 


Faint patches of NLC in clear sky seen 
from NE England, and above haze 
layer (to 5°) from Edinburgh. 


NLC visible to low latitudes over 
UK—seen first from Norfolk just 
above and parallel to N_ horizon; 
developed to herring-bone formation. 
Stacking of S-shaped bands seen 
Northolt 0140. NLC visible very 
brightly 0145-0245 from Plymouth. 
Intense glow in NNE reduced to 
bright fibrous strands as sun rose. 
Southerly stations reported orange 
tinge at N edge of cloud field. Maxi- 
mum extension around 0200. 


Display visible only in northern and 
central Ireland, possibly due to haze 
interference in many areas. When 
clear of haze, NLC forms very bright, 
e.g. around 0200 as seen Aldergrove. 


In Aberdeen two observers, inde- 
pendently, suspected presence of 
NLC visible through haze at high 
elevation. 


NLC visible Denmark with gradually 
increasing brightness. 


Prolonged, extensive and bright dis- 
play of NLC widely reported Scotland, 
N and central England, though seen 
dimly through haze in some areas. 
N-S direction of bands with fibrous 
cross structure. At 0200 bright band 
stretched almost to zenith over 
Leuchars, and at Boulmer 0220-0230 
NLC reached high elevation to east 
and west of station. Whirl formation 
reported at 0145-0215. Elevation of 
N edge unreliable due to haze. 


Band of NLC visible NW-NE from 
Edinburgh. 


NLC seen, as soon as sky darkened, 
to occupy large area of twilit hemis- 
phere—to zenith as seen from W. 
Yorks. and Berkshire—until obscured 
by moonlit tropospheric clouds. NLC 
photographed 2210 W. Yorks. 


56°N 03°W 

S5S6°N 04-5°W 
56°N 03°W 

55:5°N 07:5°W 
55:5°N 01-5°W 
55:5°N 12°E 
SS°N_ 14-5° 


55°N 
54-5° 


54°N 
53°N 


56°N s 
54:5°N 06°W 
54°N 04-5°W 
53-5°N 07-5°W 


56°N 03°W 
55-5°N 01-5°W 


54-5°N 06°W 
53°N 08°W 
52:5°N 0-5°E 
52°N_ O1°E 


51-5°N 02°W 
51-5°N 0:5°W 


S1°N_  01-5°W 
50-5°N 3-S°W 
50-5°N 05°W 
55-5°N 07-:5°W 
54-5°N 06°W 
53°N 08°W 


S7T?7N 02°W 


55°5°N 12°E 


57-5°N 03-5°W 
S57°N _ 02°W 
56°5°N 03°W 


56°N 03°W 
55-5°N 01-5°W 


53°N 05°W 


56°N 03°W 


54°N 02°W 


53°N_ _0-5°W 
51-5°N 0-5°W 


_ 


—— 
COA CINK IS 


+ 


360 
345-015 


~ 360-040 


355-060 
022-045 


350-020 


360 
340-040 
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TABLE I—continued 


Date— a Station : Max. Limiting 
night of Times Notes Position Time elev. azimuths 
UT UT degrees 


12/13 2250 NLC visible through breaks in tro- 56°N 03°W 2250 13+ 045 
pospheric cloud (to 13°). 


15/16 2330 Low elevation NLC seen after clearing 57:5°N 03-S°W 340-015 
0050-0215 of tropospheric cloud; thin delicate 57°N 02°W 340-010 
structure. S6°N 03°W 340 
55-5°N 04-5°W 340-010 
5S°N 04-5°W 


Short-lived appearance of NLC in E  57-5°N 03-5°W 080-090 
in break in tropospheric cloud. 


Compact formation NLC from NE 59°N 09°E 
horizon to 40° elevation seen Norway. 


2346-0220 Whirl formation seen from northerly 57°N 02°W 
stations; more general reporting of 
bright bands spreading to high eleva- 
tion over N Scotland and visible from 56°N 03°W 
central England. NLC photographed 56°N_ 10°E 
Aberdeen 2353-0219. Observation 
from Denmark hampered by cloud 55-5°N 01-:5°W 
and moonlight, but brighiening NLC 
bands seen in cloud breaks. 01-5°W 


2200, 0100 Possible display of NLC bands high 09°E 
in SE seen Norway; faint NLC sus- 03°W 
pected Edinburgh. 


2335-0100 In Denmark bands of NLC seen first 10°E 
through breaks in tropospheric cloud; 
wae brightness; photographed 


2135-0105 In Denmark first seen faintly in fine 
observing conditions to 12° but dis- 
appeared 2205. Visible again 2400 
with finely striated structure fading 
later into light sky. 


31 July/i Aug. 0040-0200 Very bright display of NLC with 
typical band and billow formation of 
degree of brightness and extent not 
previously seen by the observer at this 
time of year. Series of photographs. 


2/3 Aug. 2145 Possible weak display of NLC repor- 
0145-0305+ ted at earlier time from Norway. 

From 0145 to about 0300 bright NLG 
bands visible from two of most 
northerly Met. stations and Fair Isle 
lighthouse to remarkably high eleva- 
tion in conditions almost clear of 
tropospheric cloud. 


Suspected NLC 
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551.577.37(41-4) 


EVALUATING THE PROBABILITY OF HEAVY RAIN 
By M. C. JACKSON 


SUMMARY 


Some of the weaknesses of the traditional Bilham methods for estimating rainfall depth, 
duration and return period are discussed. It is shown how much improved estimates for 
anywhere in the United Kingdom can be made using the methods described in Volume II of 
the Flood Studies Report (Meteorological Office, 1975). 

The methods are based on maps of the 60-minute and 2-day rainfalls with return period 
5 years. This paper describes how to obtain estimates of rainfall for the same return period 
for other durations using these two basic maps, and also how to extend rainfall estimates up 
to return periods of 1000 years. A method of estimating areal rainfall from the already 
derived point rainfalls is also described. 


1. INTRODUCTION 


Engineers concerned with urban drainage, sewerage, management of river 
catchments, planning construction, etc. have for many decades requested 
information about rainfall amount and return period for different parts of the 
country and for various durations and sizes of area. (A rainfall with return 
period T years is one that has a probability of 1/T of being equalled or exceeded 
in any one year.) In the past such enquiries have usually been answered by 
reference to the Bilham equation and tables, which gave the engineer an estimate 
of the probability of return period of a specified amount of rain falling in a 
specified duration. The equation was first presented in British Rainfall (Bilham, 
1935). 

However, the amount of data used by Bilham was limited to only 120 station 
years; other weaknesses are that only one equation was derived to represent the 
whole of the United Kingdom, and that durations were limited to two hours or 
less. Consequently estimates of rainfall with return period greater than 5 years 
were often in considerable error. 

A modified Bilham equation (Holland, 1964) was more complex, adjusting 
the results when the intensity is greater than 1-25 inches per hour (32 milli- 
metres per hour). However, no other changes were made to improve estimates of 
rainfall values with return periods greater than 5 years, and, in fact, further 
rainfall values for durations from 2 hours to 24 hours were obtained by simply 
extrapolating the Bilham equation beyond 2 hours. 

Consequently there has been a need for many years both to collect many more 
data and to devise techniques for accurately estimating depth, duration and 
return period information for individual places in the United Kingdom. 

In 1969 a major project, the United Kingdom Flood Studies, was established 
by the Natural Environment Research Council (NERC). The meteorological 
studies were directed by A. F. Jenkinson of the Meteorological Office, and one of 
the main products of the three years’ work (1970-73) was Volume II (Meteoro- 
logical Studies) of a comprehensive Flood Studies Report (Meteorological 
Office, 1975). For this work 3000 station years of data from autographic 
recording rain-gauges and 40 000 station years of daily rainfall data were collec- 
ted and analysed. Most of the methods described by A. F. Jenkinson and M. C. 
Jackson in the Flood Studies Report, and by Jenkinson (1976), and described 
briefly in this paper, can be used to give estimates of rainfall and return period for 
different parts of the country and for various durations and various sizes of area. 
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More recently some of the calculations have been made by computer, and Keers 
and Westcott (1977) describe these computer programs. 


2. MAPPING THE MAGNITUDE OF RAINFALL EVENTS WITH RETURN PERIOD 
5 YEARS 


Key parameters were needed which could be mapped to give fundamental values 
on which to base estimates of rainfall amounts with long return periods (as well 
as being used to standardize the very large amounts of information obtained 
from the data). The key return period chosen was 5 years, and events for all 
durations with a return period of 5 years (sometimes called the M5) were linked 
together. 

A return period of 5 years was chosen as a standard since it satisfies two basic 
requirements: (i) the larger the return period which is used as a standard, the 
better will be the estimation of events with even larger return periods, (ii) the 
smaller the return period, the better the chance of estimating its value accurately. 
The return periods satisfying (i) and (ii) are between 2 and 20 years, and 5 years 
was chosen for theoretical reasons. (The magnitude of the M5 rainfall event can 
be estimated by the reader by calculating the geometric mean of the top half of 
an ordered set of annual maxima—see Flood Studies Report, Volume II, 
Chapter 2.2.) 

In order that M5 events can be estimated at any specific place for all durations, 
the magnitude of the MS event for each of two durations at that place are 
calculated, and then, with a relationship derived from all the data, the magni- 
tudes of the M5 event for all the other durations can be estimated (Section 3). 

The two key durations chosen for M5 mapping were 60 minutes (Figure 1) 
and 2 rainfall days (Figure 2), after Flood Studies Report, Volume II, Figures 
3.6* and 3.3 respectively. Sixty minutes is the shortest duration for which rainfall 
is usually tabulated One rainfall day (0900-0900 GMT) would seem the most 
logical duration for the second key event, with many cyclonic disturbances 
giving heavy rain which lasts for about a day. However, the rainfall event is often 
cut in two by an 0900 observation: this problem is overcome by choosing a 
duration of 2 rainfall days. The 2 rainfall-day event is smaller than the 48-hour 
event because a 48-hour event can commence at any hour of the day and 
continue for 48 hours, while the 2 rainfall-day event, although lasting for 48 
hours, can only start at the 0900 morning observation. The ratio of the 48-hour 
event to the 2 rainfall-day event with the same return period is approximately 
1-06—see Flood Studies Report, Volume II, Chapter 3.3.3. The same happens 
with rainfall data for 60 minutes compared with a clock hour and with data for 
30 days compared with a calendar month. 


3. ESTIMATION OF THE MAGNITUDES OF RAINFALL EVENTS WITH RETURN 

PERIOD 5 YEARS FOR ALL DURATIONS AT ANY PLACE 
Autographic rainfall data from recording rain-gauges in the United Kingdom 
were analysed extensively to give the table below. Given the values of the M5 
rainfall for duration 60 minutes and 2 rainfall days, and the ratio of the two 
values, Table I allows estimates to be made for durations from 1 minute to 
48 hours. Table I is based on Flood Studies Report, Volume II, Tables 3.6 and 
: A 





a La map originally published as 3.6 was inserted in error, and a corrigendum has been 
issued. 
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FIGURE 1—MaAP OF RAINFALL AMOUNT FALLING IN 60 MINUTES WITH 
RETURN PERIOD 5 YEARS (M5) 
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FIGURE 2—MaAP OF RAINFALL AMOUNT FALLING IN 2 RAINFALL DAYS WITH 





RETURN PERIOD 5 YEARS (M5) 
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TABLE I—METHOD OF ESTIMATING M5 RAINFALL FOR ALL DURATIONS 


Duration (d) in minutes Duration (D) in hours 
5 1 15 30 120 2 4 6 12 24 48 


“4 120 92 106 
"3 123 89 106 
ts 126 86 106 
‘Z 130 3 83 106 
“2 134 80 106 
“1 139 77 =106 
pz 149 72 


r = ratio of M5 rainfall for 60 minutes to MS rainfall for 2 days. 
x = ratio (per cent) of M5 rainfall for d minutes to MS rainfall for 60 minutes. 
X = ratio (per cent) of MS rainfall for D hours to MS rainfall for 2 days. 


The relationship of M5 rainfall with duration at any one station fits a simple 
numerical model, described in Flood Studies Report, Volume II, Chapter 3.6: 


In J—In In, = —n In (1 + BD) 


where J is the rainfall intensity (mm/h), 
I is the rainfall intensity for very short durations (viz. 15s), 
D is the duration in hours, 
B and n are factors which are constant for any particular location. 


Values of J, range from 175 mm/h at stations in south-east England to less 
than 100 mm/h in Scotland, while n ranges from 0-77 in south-east England to 
less than 0-60 in mountainous parts of Scotland, and B ranges from 15 in the 
drier parts of the country to more than 30 in wetter, more mountainous regions. 


4. ESTIMATION OF THE MAGNITUDE OF EVENTS WITH OTHER RETURN PERIODS 


It is usually necessary to know the magnitudes of rainfall events with return 
periods longer than 5 years. The rate of increase of rainfall with return period 
varies from station to station, reflecting the details of the heaviest falls at each 
station. However, when dozens of stations are gathered together into regions, 
the regional average rate of increase becomes more or less the same, wherever 
the region is. Because of this, rainfall data for the whole country were used to 
establish the ratios of the 7-year event to the 5-year event, M7/MS, called the 
growth factor. 

In Figure 3, drawn from Table 2.7 of Flood Studies Report, Volume II, the 
growth factor is shown against the magnitude of the MS rainfall event: a single 
set of curves is sufficient to estimate rainfalls with return periods other than 5 
years in widely different places, such as Cambridge and Snowdonia, even though 
their other rainfall characteristics are so different. For example, if the 6-hour 
MS event in Snowdonia and the 24-hour MS event in Cambridge are of the same 
magnitude, then the events with longer return periods, but still of these same 
durations, are also of the same magnitude. 

In fact there are some small differences in growth factor from region to region; 
the differences between north and south are just large enough to justify treatment 
as two separate regions; England and Wales (Figure 3), and Scotland and 
Northern Ireland (not given)—see Flood Studies Report, Volume II, Chapter 2.3. 
Growth factors were calculated for return periods of up to 1000 years, with 
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FIGURE 3—RATIO OF 7-YEAR RAINFALL EVENTS TO 5-YEAR EVENTS PLOTTED 
AGAINST 5-YEAR RAINFALL (M5) AND RETURN PERIOD; SUITABLE FOR ANY 
POINT IN ENGLAND AND WALES 


rough estimates for as much as 10 000 year::: quite clearly the larger the return 
period, the less confidence we have in the esiimate. For further discussion of 
this, the reader is referred to Flood Studies Report, Volume II, Chapter 2.3. 


5. MAGNITUDE OF A RAINFALL EVENT OVER AN AREA 
Engineers are often more interested in the magnitude of an areal rainfall event 
over a given catchment than in a point event such as those just described. To 
meet this requirement, areal reduction factors have been calculated, which 
depend on the size of the area and the duration of the event. When an areal 
reduction factor is multiplied by the value for the point rainfall event, the 
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product gives the magnitude of the areal rainfall event with the same duration 
and the same return period. It must be emphasized that this type of areal reduc- 
tion factor is used both here and in the Flood Studies Report; it is not applicable 
to individual point-centred storms, but is for catchments which are independent 
of storm location. 

The factor is smallest for small durations and over large areas, approaching 
unity as the area becomes smaller and the duration larger. Some point-centred 
reduction factors for the country were produced by Holland (1964), but the 
factors needed here for this kind of calculation are very different from Holland’s. 
A graph for areas of size between 10 and 2000 km? is shown in Figure 4, which is 
based on Flood Studies Report, Volume II, Chapter 5.2. 
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FIGURE 4—AREAL REDUCTION FACTORS: RATIO OF THE AREAL RAINFALL TO 
THE POINT RAINFALL OF THE SAME DURATION AND RETURN PERIOD 


6. SIMPLE EXAMPLES 

Example 1 

An engineer asks for the magnitude of the 7-hour, point rainfall events in 
Wellingborough, Northants with return periods of 5, 10, 20 and 100 years. 

Wellingborough is found on a map and the national grid reference noted 
(NGR SP 8968). The 60-minute M5 map is examined (preferably a larger version 
than Figure 1—copies are available in the Meteorological Office) and a value of 
about 19-2 mm read off at the appropriate grid reference. Similarly the 2-day 
map is examined (preferably a larger version than Figure 2—copies are available 
in the Meteorological Office) and the value of about 44-8 mm read off. Using the 
ratio of 60-minute to 2-day rainfall, 19-2/44-8 = 43 per cent, a value of 67 per 
cent can be interpolated from Table I for the ratio of 7-hour to 2-day rainfall. 
This gives a value for the 7-hour M5 as about 30-0 mm. 
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Further estimates can now be made of the magnitude of 7-hour rainfall events 
with longer return period. Figure 3 is a general diagram for England and Wales, 
showing the ratio of the MT event to the MS event, plotted against M5 values. 
Values are read off for M5 = 30-0; M7/MS is 1:21 when 7 = 10, 1-41 when 
T = 20 and 1-97 when T = 100 years. Multiplying the M5 value of 30-0 mm by 
these factors gives 36-3, 42-3 and 59-0 mm for the magnitudes of the 10-year, 
20-year and 100-year 7-hour rainfalls respectively at Wellingborough. 


Example 2 


An engineer asks for the magnitude of the 30-minute, 60-minute and 120- 
minute rainfall events with return period of 50 years over a catchment of area 
15 km? centred near Axminster, Devon (NGR SY 2998). 

Axminster is found on the maps of 60-minute M5 and 2-day MS as in Example 
1, the values of 19-0 mm and 63-0 mm are read off. Using the ratio of 60-minute 
to 2-day rainfall, 19-0/63-0 = 0-30, a value of 78 per cent can be interpolated 
from Table I for the ratio of 30-minute to 60-minute rainfall, and a value of 38 
per cent for the ratio of 120-minute rainfall to 2-day rainfall. These give the size 
of the point rainfall events for 30-minute, 60-minute and 120-minute events as 
14-8, 19-0 and 24:0 mm. 

The ratio of M50 to M5 for M5 values of 14-8, 19-0 and 24-0 mm are found 
from Figure 3 to be 1-70, 1-725 and 1-725, which gives the M50 events as 
14-8 x 1-70 = 25-2, 19-0 x 1-725 = 32:8 and 24-0 x 1-725 = 41-4 mm for 30, 
60 and 120 minutes respectively. 

Reference to Figure 4 shows the areal reduction factors for 15 km? to be 0-86. 
0-89 and 0-92 for 30, 60 and 120 minutes. Multiplication by 25-2, 32:8 and 
41-4 mm gives the areal falls with return period of 50 years as 21-7, 29-2 and 
38-1 mm in 30, 60 and 120 minutes respectively.* 
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* Note that the application of the methods to larger catchment areas, especially with 
~_— —— requires more care. Examples are given in Flood Studies Report, Volume 
, Chapter 8. 
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The measurement of airborne particles, by R. D. Cadle. 150 mm x 230 mm, 
pp. xi + 342, illus. John Wiley and Sons Ltd, Baffins Lane, Chichester, 
Sussex, 1976. Price £13-85. 


This book is written for scientists and engineers who have to sample airborne 
particles and to measure their concentration, size and shape. It gives illustra- 
tions and diagrams of instruments, the basic theory underlying them, and 
practical advice on their use. The emphasis is on commercially developed and 
marketed instruments, and the book will be helpful to anyone lucky enough to 
have dollars to spend on the generally high-class American instruments in the 
air pollution field. The treatment is clear and, although sometimes abbreviated, 
is adequate for the book’s purpose and is backed by the standard references. 
There are no obvious omissions. Some readers may consider that space devoted 
to manufacturers’ photographs and to the basic theory of such well-known 
instruments as the optical microscope could have been put to better use. 

The book is not intended as a textbook on aerosols. Meteorological applica- 
tions are only briefly referred to. There is a short discussion on the characteris- 
tics and evolution of the size distribution of the atmospheric aerosol but there 
is no discussion on the condensation of water or ice on particles. Holography 
receives 21 lines of text and two references. Two pages, inserted in a chapter on 
Optical Microscopy, are devoted to the use of particles to trace atmospheric 
motions. 

Anyone starting work with aerosols or pollution will find this a useful intro- 
ductory handbook to methods, instruments and to the literature, but meteorolo- 
gists will need to look elsewhere for full discussion of physical principles and 
for meteorological applications. 


A. C. CHAMBERLAIN 


Physical principles of micrometeorological measurements, Developments in 
Atmospheric Science, 6, by P. Schwerdtfeger. 240 mm x 170 mm, pp. ix + 
113, illus. Elsevier Scientific Publishing Co, PO Box 330, Amsterdam, The 
Netherlands. Price: Dfi 75. 


There are a number of standard textbooks on micrometeorology but most 
provide little information on the use of instruments which make the specialized 
measurements necessary in micrometeorological investigations. Professor 
Schwerdtfeger’s slender volume attempts to fill this gap: in his own words, the 
book ‘has been designed to emphasize the physical basis of precise meteoro- 
logical measurements, especially those necessary in determining processes 
relevant to the atmosphere’s behaviour close to ground level’. The material of 
the book is based on a university course and a substantial part is occupied by a 
description of about thirty experiments intended to be carried out by students. 
The experiments themselves are used both to illustrate the practical aspects of 
measurements and as a means of developing some of the principles involved. 

Chapter 1 is concerned with air temperature measurements and the transfer 
of sensible heat and deals, inter alia, with topics such as thermal inertia and 
heating by radiation of thermometers, and heat transfer from plane surfaces. 
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Chapter 2 discusses the determination of solar and terrestrial radiation and 
related topics including estimation of the solar constant, and albedo. The book 
continues with a brief discussion of atmospheric pressure, followed by a section 
on the measurement of humidity. Chapter 4 introduces the concept of eddy 
transfer, following on from wind speed measurements and interpretation of the 
vertical gradient of wind speed. The fifth chapter deals with heat conduction in 
the ground. The final chapter, which is on the modelling of thermal processes by 
electrical analogues, reflects a particular interest of the author but seems to be an 
inappropriate contribution in such a short book. 

The general level of presentation is elementary enough for the pace of the 
volume to be describable as ‘leisurely’, but this has resulted in a rather super- 
ficial treatment of many topics, and significant omissions. As a result, the book 
fails to live up to the assertion on its jacket that it is of great value as a textbook 
and reference work for meteorologists. On the other hand, if its remarkable 
price can be ignored, it provides a useful introduction for newcomers to the 
field. However, potential readers should be warned that c.g.s. units are used in 
many places in the text. Also, the sections on wind profiles and turbulent 
transfers are sprinkled with loose statements which can scarcely provide the 
sound foundation of knowledge required by the student reader: for example, the 
implication that only in neutral conditions is the mean vertical velocity near the 
ground essentially zero (pages 59 and 60), that friction velocity is the geometric 
mean of the horizontal and vertical velocity fluctuations (implying perfect 
correlation between them) (page 60), that atmospheric turbulence is isotropic 


(page 60), that the wind velocity vanishes at a height above the ground equal to 
the roughness length, and that the roughness length over developed waves is 
dependent on wind speed (for which there is little evidence) (page 61), and that an 
array of only 100 buckets arranged on the ground upwind of a profile mast will 
produce, in neutral conditions, a fully modified logarithmic profile of wind speed 
up to a height of at least 4m. Blemishes such as these significantly reduce the 
value of what is otherwise a useful addition to the micrometeorological literature. 


N. THOMPSON 


NOTES AND NEWS 


Retirement of Mr R. Murray 


Mr Robert (Roy) Murray, who retired on 28 April 1977, joined the Office as a 
Technical Officer in September 1939 after graduating at Edinburgh University 
with an M.A. degree in Mathematics and Natural Philosophy. After a short 
initial training period in London he went to Headquarters No. 221 Group 
(Coastal Command), Royal Air Force Donibristle, later Pitreavie, for forecasting 
duties. Early in 1941 he was posted to the Isle of Islay and remained there for 
about a year before proceeding to Northern Ireland for relief duties at a number 
of RAF Coastal Command units. His service in Northern Ireland concluded 
with a spell at RAF Headquarters and in March 1943 he was commissioned in 
the RAFVR and sent to India, joining HQ No. 221 Group initially in Bengal. 
Fit Lt Murray proceeded with this Group to Imphal, Manipur, in November 
of that year. He remained at Imphal for a year or so and throughout the 
Japanese siege. He then became the Senior Meteorological Officer with the 
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mobile combined Headquarters 14th Army and 221 Group and travelled by 
stages through Burma to arrive in Rangoon in mid 1945. On the termination of 
hostilities he filled the post of Senior Meteorological Officer Burma and was 
responsible for a number of units scattered throughout Burma, Siam and Indo- 
China. He was mentioned in dispatches at this time. In February 1946 Sqn Ldr 
Murray (as he then was) was repatriated and demobilized, and took up civilian 
duties in the Central Forecasting Branch (Met O 2) at Dunstable. He became a 
Senior Scientific Officer in September of that year and in 1950 he moved to the 
Forecasting Research Branch (Met O 21) where he did some first-class work in 
pioneering the understanding of the dynamics of jet streams which was reported 
in a classic paper. Promotion to Principal Scientific Officer followed in 1953 
with a return to operational forecasting duties. In August 1956 he was posted to 
Aden to become the Chief Meteorological Officer, MEAF. Two years later he 
returned for a brief spell to a research post in upper-air climatology at Harrow 
(Met O 13) and subsequently went overseas again in mid 1959 to Cyprus to 
occupy the Chief Meteorological Officer post at Headquarters NEAF where he 
remained for the next four years before returning in August 1963 to the Synoptic 
Climatology Branch located at Bracknell. He spent the next decade in that 
Branch contributing to research in forecasting for periods of one month to a 
season. It is interesting to record that the first 30 day forecast issued to the 
public (for December 1963) correctly predicted a major change from an ex- 
ceptionally warm November to an exceptionally cold December and that this 
was largely due to his confident prediction of an imminent block. Throughout 
his career in extended-range forecasting his approach was essentially a practical 
one based on a desire to establish the technique on a sound and objective basis. 
During this time he published some 30 papers on the detailed climatology of the 
British Isles, developing objective indices and rules for forecasting monthly and 
seasonal temperature and rainfall. He was awarded the L. G. Groves Memorial 
Prize for Meteorology in 1970 (jointly with R. A. S. Ratcliffe) for work relating 
ocean temperature to atmospheric circulation anomalies. 

In May 1973 Mr Murray was promoted to Senior Principal Scientific Officer 
and took charge of agricultural meteorology and hydrometeorological work as 
an Assistant Director. He remained in this post until his retirement. During 
this period he successfully led efforts to establish and introduce computer-based 
techniques which enabled the services provided to the community in general, 
and the water and agricultural industries in particular, to be improved in 
quality, quantity and range. In the last year of his service he was actively 
involved in many aspects of work resulting from the 1975/76 drought. 

There can be few left in the Service who have had as varied a career. It was 
not until the last 18 months that I came into day-to-day contact with him, and 
I would like to record my personal thanks for his informed advice and willingly 
tendered help. We wish Roy and his wife a long and active retirement. 


N. E. RIDER 
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